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|^ I. INTRODUCTION 

o 

Novel symmetry types of the energy gap function A(k, w) in spin-singlet superconductors were discussed recently 1 in 
the context of high- T c superconductivity, 2 ' 3 in order to investigate ways of constructing a superconducting state which 
is not inhibited by the the presence of strong short-range repulsion. An anisotropic superconducting gap parameter 
"j-^ ( "p-wave" or "d-wave" ) which changes sign as k moves around the Fermi surface is the simplest way of achieving such 
a state. An alternative posssibility was proposed some time ago by Berezinskii, 4 who pointed out in the context of 
3 He that it is permissable to have triplet (5 = 1) pairing in which the gap function A(k, w) is even in momentum 
l k (for example "s-wave"), provided it is odd in frequency w. Recently, this suggestion was extended to the singlet 
1 O case, 1 where it was argued that a new class of singlet superconductors can exist for which the gap function is an odd 
£h function of both k and w. Such a superconducting state is thus odd under both parity and time reversal and, as in 
the case discussed in Refs. 2, 3, the effect of short-range repulsion is suppressed, here because the equal-time pair 
amplitude is zero. 

An odd-in-frequency solution to the Eliashberg equations can only exist for certain forms of the pairing kernel 
which scatters a pair from k, — k; w, — w to k', — k'; w', — w'. In the singlet case, to obtain an odd w gap parameter, 
the kernel must have a contribution which is odd in both k and w. The problem of how an appropriate interaction 
might arise has been discussed in a recent paper 5 where a specific example based on spin-fluctuation exchange in 
the two-dimensional Hubbard model was discussed. If an interaction admits an odd solution it still must be checked 
whether it only describes a metastable phase with free energy higher than that of an even-frequency solution. 

In the present paper, we explore another issue which arises because there is an instability which was not mentioned 6 
in all the known translationally invariant odd-gap solutions considered in Refs. (1,4). Namely, because of the symmetry 
properties of the odd-gap pairing amplitudes, there is a change in sign of the anomalous parts of BCS coherence factors. 
A consequence is that a calculation of the Meissner effect in lowest order 7 gives an opposite sign to that of the usual 
BCS case so that the superfiuid density appears to be negative. Thus, although the free energy of the uniform odd 
gap state is less than that of the normal state, 1 the system prefers a non-uniform distribution of supercurrents, i.e. 
some sort of vortex state which could in principle be determined in a calculation which goes beyond the linearized 
Eliashberg equation which is used to find T c . 

In this paper, we shall show how a stable Meissner state may be achieved by the introduction of a composite 
condensate which is characterized by the joint condensation of a Cooper pair and a density fluctuation. Our point 
of view is motivated by the question: If the anomalous Green's function vanishes at equal times, as it does for an 
odd gap, what is the appropriate macroscopic thermodynamic order parameter? For example, in the usual BCS case 
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it is the expectation value F(r,t;r',t' — > t) = (ip(r, t)ip(r' , t)) for |r — r'| < £, where £ is the BCS coherence length. 
Since this vanishes in the odd case (it is odd in t — t'), we are led naturally to consider the quantity [dF(r, t; r', t')/dt] 
for t — > t' as the equal-time correlation function to be used as order parameter. Given a pairing hamiltonian, this 
quantity can be obtained from the equations of motion and it consists of the expectation value of the product of a pair 
operator and, depending on the nature of the pairing interaction, a spin or charge density fluctuation. This quantity 
is of the form 

<M(r,*).V(M)V(M'^)}, (1.1) 

with an appropriate spin structure which will be discussed below. The operator M contains a spin or charge density 
fluctuation: 

M(r,t) = J dTdr'^(r',t + T)iP(r',t + T) £>(r - r', t), (1.2) 

where V is a generalized non-local retarded interaction. 

In Section II, we present a symmetry analysis of the Green's functions which leads to this new type of superconductor. 
We show that for the spin-singlet (triplet) case, the only symmetry requirement is that the gap function A(k, w) be 
even (odd) under simultaneous change of sign of momentum and frequency. Consequently, apart from the conventional 
even-frequency gaps which are even (odd) in momentum for singlet (triplet) pairing, there can be solutions of the gap 
equation which are odd (even) in momentum and odd in frequency in the singlet (triplet) case. 

In Section III, we give an illustration of how a composite operator arises naturally as the order parameter for an 
example system. For purposes of calculation, we introduce a model hamiltonian which can lead to odd-frequency 
pairing and we derive the gap equation near the critical temperature in terms of the parameters of the model. We 
find that the strength of the interaction must exceed a critical value in order that the odd solution exist. This is a 
rather general property of odd-gap solutions 4 ' 3 and is in contrast to the BCS case where the logarithmic behavior of 
the gap equation insures a non-trivial solution no matter how weak the coupling. In our case, the kernel of the gap 
equation, being odd in w, must vanish at zero frequency if the kernel varies smoothly in this vicinity; this removes 
the logarithmic divergence. 

In Section IV, we calculate the Green's functions in the composite operator condensate and show explicitly that 
the superfiuid stiffness is positive. Section V contains a discussion of the Josephson effect between odd-gap and even- 
gap superconductors. In Section VI we calculate the NMR relaxation rate and in Section VII the effect of impurity 
scattering for the odd-gap composite superconductor 

Subsequent to the singlet odd-gap discussion of Balatsky and Abrahams, 1 several authors, in particular Kivelson and 
Emery, 8 Balatsky and Bonca, 9 and Coleman, Miranda and Tsvelik, 10 have proposed other scenarios for odd pairing 
involving composite operators. These examples indicate that the odd-frequency-gap superconducting correlations are 
relevant for some physical models. Although a robust realization of odd-gap superconductivity remains to be found, 
in what follows we consider a simple model in which some properties of the odd gap superconductors, such as the 
Meissner effect, NMR and Josephson effects, can be analyzed. 



II. SYMMETRY OF THE GAP 



We introduce the anomalous Green's function 

T(k,u n ) = V / ' dre-^- T (T rCkia (r/2) C _ ki/3 (-T/2))^ a/3 (k), (2.1) 

where u>„ = (2n + l)-7rT and <^ a| g(k) determines the spin symmetry of the pairing amplitude and hence the gap 
function. For example, for singlet pairing, <f> a p = i((j y ) a p/2, where <r y is the Pauli spin matrix. For triplet pairing 
(as in 3 He, for example) <j> a p(k) = {i(T y <j) a p ■ d(k), where d(k), a vector in spin space, is a linear combination of the 
three Y£ ra (k) for I = 1, each of which corresponds to one of the m z = ±1, states. For example, the form d^ ± id,, 
gives the m = ±1 component of the triplet. Similar forms hold for the anomalous self energy W(k, w„) and the gap 
function A(k, w„). n 

The only constraint on the possible symmetry of ^(k, u>„) and VF(k, u>„) follows from the anticommutativity of the 
ij> operators in T. Using this and the antisymmetry of the Pauli matrices we immediately get, for the singlet case: 11 

T(k,u n ) = ^(-k,-w„), A(k,w„) = A(-k,-w„) (2.2) 
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and for the triplet 4 

T(k,u n ) = -^(-k,-w„), A(k,w„) = -A(-k,-w„) (2.3) 

From these relations, it can be seen that the conventional even spatial parity singlet (conventional metals) and odd 
parity triplet ( 3 He) states have gap functions which are even in Matsubara frequency. Then the equal time anomalous 
Green's function T is nonzero and this leads to the usual off-diagonal long-range order, ODLRO. However it has been 
shown 1 ' 4 ' 5 that gap functions which are odd in Matsubara frequency are admissible. In the odd case, Eq. (2.2) is 
satisfied with singlet superconducting pairing with a pair amplitude which is odd in both k and w„: 

T{k,u n ) = -T(-k,u n ) = -T{k,-u n ) 

A(k,w„) = -A(-k,w„) = -A(k,-w„). (2.4) 

Similarly, for the triplet case, Eq. (2.3) is satisfied for triplet superconducting pairing with a gap which is even in k 
and odd in w„: 

T{k,u n ) = T(-k,u n ) = -T{k,-u n ) 

A(k,w„) = A(-k,w„) = -A(k,-w„). (2.5) 

Eqs. (2.4, 2.5) show that the spin-singlet (triplet) gap is described in terms of an odd (even) orbital function, while, 
at the same time, the spin function is odd (even). There is no violation of the Pauli principle because the equal-time 
gap function, or pair amplitude, for singlet or triplet vanishes since the gap is odd in oj n . 12 



III. EQUATIONS FOR COMPOSITE ODD SUPERCONDUCTIVITY 

A. Model for Composite Operator 

For simplicity, we consider the m = 1 triplet Berezinskii 4 state. The conclusions are similar for the singlet case, 
but the notation is somewhat more cumbersome. From Eq. (2.1), we write the anomalous Green's function as 

T{k, t) = ^(T r c kja (r) c_ kj/3 (0))(i^) a/3 • d). (3.1) 

a/3 

which, for d = — (d^ + id y )/2 (for example), reduces to 

T(k,r) = (T r c k , T (r) c_ k , T (0)}. (3.2) 

We are considering the odd gap case in which T{t) is odd in r and even in k, an even parity triplet odd gap. Since 
T(t) is odd, it vanishes at r = 0. However, its T-derivative is even and is generally non-zero at r = 0. 13 Therefore, 
we study the anomalous amplitude 

^=^^ = (T,^ C _ kT (0)). (3.3) 

The behavior of T is deduced from the equations of motion. The hamiltonian is 

n = H c + H s + H x , (3.4) 

where H c is the kinetic energy of the conduction electrons, Hs is the hamiltonian for some spin or charge excitation 
and Hx couples the electrons to the excitations. If one integrates out the excitations to produce an effective retarded 
(perhaps spin-dependent) electron-electron interaction as in the Eliashberg 14 treatment of superconductivity, then T 
can be evaluated by commuting c^(t) with the effective Hamiltonian. In that case, T has a term of the form of Eq. 
(l-l)- 

In order to motivate the structure of the odd composite condensate, we consider a spin fermion model in which the 
conduction electrons are coupled via Hx to localised spin modes which arise from a second band of electrons. Thus, 
Hx has the form 
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(3.5) 



where S 8 - is a spin 1 excitation at site i of the lattice. In what follows, we shall assume that both the possibility of 
Kondo screening (for J > 0) and the induced ("RKKY") interaction between the local spins can be neglected in the 
present context. 

As in Eq. (3.3), the pairing is studied by means of the anomalous amplitude Ti{r) which is determined by c 8 - = 
[H c + H x ,a]: 

Mr) = - X)<.-i(2Vc iT (r)c,- T (0)) - Ja w • L </3 (t), (3.6) 

3 

where t,j is the lattice kinetic energy matrix element and L ! - | g(r) is the composite anomalous amplitude which 
characterizes the condensate. L is a vector in spin-space and for the m = 1 odd condensate, it has the form 
Lj-^(t) = (T r C8 | g(r)c8i-(0)Sj(r)). We shall be interested in a spin-one condensate in which a spin-zero Cooper pair is 
coupled to a spin-one excitation (magnon) S 1 * 1 . Therefore, we shall study 

L(l, 2, 3) = L a/3 (1, 2, 3)e a/3 /V2 = (T r c la (T 1 )c 2/3 (T 2 )S 3 (T3))£ a/3 /V2, (3.7) 
where 1,2,3 label space and imaginary time points and e a fi is ic^a. 



B. Reduced Hamiltonian and Mean-Field Equations 

Our aim is to calculate properties of the odd pairing state which depend on the structure and symmetry of the 
composite order parameter. For this purpose, we start with the simplest possible weak-coupling model which gives rise 
to a condensate characterized by a non-zero value of L, the expectation value of the composite operator discussed in 
the previous subsection and which we can solve analytically. In direct analogy to BCS theory, 15 we write a "reduced" 
hamiltonian 

H re d = H + y^cjq C !<3 S » • Vij CjyCj 6 S^ (CapCty). (3.8) 
ij 

Here V,j is an attractive (i.e. V < 0) short-range, instantaneous interaction which mediates the condensation and Ho 
contains the kinetic energy of quasiparticles and the non-interacting hamiltonian for the spins. 

We should stress that the above hamiltonian is only an instantaneous approximation to a dynamical six-point 
interaction. A dynamical theory which takes into account the retardation in L(l, 2, 3) can be constructed, for example 
starting with the interaction Hx of the previous subsection [Eq. (3.5)] 16 ' 17 and perhaps other couplings as well. 

For the reasons stated above, we work with H re d. Such a reduced hamiltonian might be obtained as an approxi- 
mation to the repeated interparticle scattering of two quasiparticles with a spin, i.e. the Faddeev scattering of two 
quasiparticles with a spin fluctuation. Alternatively, if we take the simple interaction of the previous subsection and 
add an electron-electron interaction {7, then the irreducible six-point coupling of H re 4 could arise as shown in the 
diagram of Fig. 1. 




FIG. 1. Six point coupling arising from electron-electron interaction (U) and electron-localised spin coupling (J). 
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FIG. 3. Linearized gap equation for composite condensate 



The gap equation differs from the BCS case by the presence of the spin propagator in the integral. It is essential to 
recognize that the boson described by V is not what is responsible for forming the Cooper pair part of the condensate. 
The interaction for the condensate here is the instantaneous V . We shall show in the next subsection that for a simple 
choice for the propagator T>, provided V exceeds a minimum coupling, there is an instability which determines T c . 
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C. Behavior Near the Transition Temperature 



In contrast to the BCS case, the mean-field hamiltonian, Eq. (3.9), is not quadratic in the fields so that an immediate 
diagonalization is not possible. Instead, we shall assume that the transition is second order everywhere so that we may 
construct a perturbation theory in A, which should be valid both near T c and near the minimum coupling, where A is 
small. This approach is similar to that used by Gor'kov in the microscopic derivation 7 of Ginzburg-Landau equations. 
In this way, we can obtain the transition temperature and various response functions close to the transition. The 
elementary vertices A, A are shown in Fig. 4. 





FIG. 4. Elementary anomalous vertices for perturbative treatment. 

We begin with the anomalous condensate amplitude L 8 (t). To third order, it is given by the diagrams in Fig. 5 
and the analytic expression is 




+ 




FIG. 5. Perturbation expansion for L. 



L(l,2,3) = J d4 g(l,4)A(4)g(2, 4)27(3, 4) 

+ J d4d5d6 g(l, 4)A(4) • Q(5, 4)A(5)£(5, 6)A(6)g(2, 6)27(4, 6)27(3, 5). 



(3.12) 



The propagators Q and 27 are for quasiparticles and spins in the normal state. The transition temperature is found 
from Eq. (3.10) and the linearized equation for L [compare Eq. (3.11)]: 



A(i) = yL(i, i, i) = v J d2 g(i, 2)A(2)g(i, 2)27(1, 2). 

Here we have taken the interaction V to be of zero range. The T c equation is therefore given by 



k,q u>,v 



(3.13) 



(3.14) 



where w, v are fermion (odd) and boson (even) Matsubara frequencies respectively. 

The precise result for T c depends on the form of the magnon propagator. It is obvious that 27(1,2) = const 
leads to an ordinary BCS expression for T c . To introduce the new physics of the composite condensate we choose 
27(q, v) = Y l[(iv) 2 — T 2 ], which could represent the dynamics of non-interacting localized Kondo-like spin resonances. 
Since 27 is independent of momentum, we carry out the momentum sums in Eq. (3.14) and find 



1 = -2tt 2 N$VT 2 P>(w + w') - 2>(w - w') 



(3.15) 
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where No is the one-spin quasiparticle density of states at the Fermi surface. We shall see below that the sums 
are dominated by their low-frequency behavior. Therefore, for simplicity, we approximate the square bracket of Eq. 
(3.15), i.e. the odd part of T>, by the separable form 



2D. 



odd 



4r- 



(3.16) 



Eq. (3.15) can now be evaluated explicitly for T c . In terms of digamma functions, the result is 

„ r , ,1 w c + iT. , .1 «r 2 



(3.17) 



where g = 2NqVY and w c is the ultraviolet cutoff of the interaction. In Fig. 6, we show a plot of T c against the 
effective coupling g. As expected, there is a minimum value of g to achieve a transition and the curve suggests that for 
a range of g above the minimum value, there is reentrant behavior. This may be checked by evaluating the right-hand 
side (RHS) of Eq. (3.14) near T c . The condensed state occurs for those temperatures, larger or smaller than T c for 
which the temperature derivative of the RHS is negative. In this way, we find that in the (smaller coupling) region of 
Fig. 6 where T c has two values T c+ , T c _ the condensate exists between them and the system is normal for T < T c _ 
or T > T c _|_. In the larger coupling region, only T c+ exists and the system is in a condensed phase down to T = 0. 
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FIG. 6. Phase diagram for composite odd gap state, showing reentrant behavior. 
Another possibility is a spatially uniform magnon propagator of the form 

2?(q, v) = - 



(ivf 



(3.18) 



In this case, the momentum delta function removes one integration and a infrared logarithmic integral (as in BCS) is 
recovered. The T c equation becomes 



2N VT? 



L (w-w') 2 + r2w + w 



w>o, lo'>o 

The result of the frequency sums gives the transition temperature as 

T c = 1.56w c e-' r / AroV . 



(3.19) 



(3.20) 



In order to find the quasiparticle spectrum in the condensed state, we examine the single-particle Green's function. 
For this, we calculate the irreducible self-energy to second order in A. The diagram is given in Fig. 7 and the expression 
is 
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E(l, 2) = (4/3)A(l) • A(2)g(2, 1)27(1, 2). 



(3.21) 




FIG. 7. Diagram for self energy. 

After Fourier transform, Eq. (3.18) takes the form 

E(k,w) = (87ri/3)7V |A| 2 T ^ V odd (u + w') 



w'>0 



With the use of the separable propagator of Eq. (3.16), we find 



iu)Z(iu)) — e ' 



(3.22) 



(3.23) 



where Q° is the normal state Green's function, e is the normal state quasiparticle kinetic energy measured from the 
Fermi energy, and 



Z(iw) = 1 - W 2 /[(iw) 2 - T 2 ], W 2 = (47r/3)|A| 2 ^r/2y. 



(3.24) 



From Eq. (3.23), we find that the density of states in the superconducting state is unchanged from that of the 
normal state, a situation which always obtains when the self-energy is momentum independent (as in the case e.g. of 
electron-phonon interactions). However, the interaction produces mass renormalization. Here, there is no gap, but 
the quasiparticle spectrum is modified to the form 



£(k)~e(k)[l- 



W 2 



e(k) 2 - r 2J 



(3.25) 



which is valid near the transition where |A| 2 and hence W 2 is small. For small excitation energy, there is a dynamical 
mass enhancement which does not appear in the density of states, but which is important in response functions. 



IV. MEISSNER EFFECT 

The Meissner effect occurs when the paramagnetic electrodynamic response is less than the diamagnetic part. The 
dc response is given by 



ji(q) = -<2ij(q) 4/(q) 

Q 8i (q) =6 ij Ne 2 /m+Q? i (q L ) 



The paramagnetic kernel is 15 

Qfi(q) = -(e 2 /4m 2 )£5>k;. f dr(T4(k + ,T)c a (k_,T)4(k'_,0)^(k' + ,0)). 

af) kk' 



(4.1) 



(4.2) 



where k± = k ± q/2. 

As in the previous section, we evaluate Q p near T c by perturbation in the order parameter A. The relevant diagrams 
are shown in Fig. 8. 




FIG. 8. Current-current correlation function for the Meissner effect. 



The analytic expression for q — > is 

<#;-G"; = 5"l A l 2 ^ k,k > ^ 2 ( k '^ 2 ( k '^')-2a 3 (k, w )g(k', w ')P(k + k', w + w '). (4.3) 

uu', kk' 

A result Q p — Q n > indicates a Meissner effect and a positive superfluid density. 

If the magnon propagator is momentum independent, there is no contribution to the Meissner kernel since the 
momentum summands are odd. Therefore, we generalize the magnon propagator to include momentum dependence. 
An extreme possibility is the case of a static, spatially uniform magnon of the factorized form 

X>(q,i/) = -tf q tf„/T. (4.4) 

When this is used in Eq. (4.3), the result is Q p — Q n > and is numerically precisely that which is found in ordinary 
BCS theory near T c . This is no surprise; the same uniform V gives the BCS expression when used in the T c equation, 
Eq. (3.14). If we spread out the delta functions, as would occur for a more realistic magnon propagator, the sign of 
Q p — Q n will not change, the superfluid density will be positive with a value intermediate between zero and the BCS 
number. 

As an illustration, we consider a magnon propagator describing a resonance, as before, but one which is spatially 
uniform: 

P(q,^) = -^ q {iv) 2_ V 2 - ( 4 - 5 ) 

When this is used to evaluate Eq. (4.3), we find a result which may be smaller than the BCS value, and, we emphasize, 
of the same sign: 

= 19.4^. (4.6) 

Qbcs r 

This result has been obtained as the leading term in an expansion in T c /T which has been assumed small. Therefore, 
there is a Meissner effect for the odd frequency composite condensate with a reduction of the superfluid density from 
the BCS value by a factor of order of the ratio of the transition temperature to the resonance width. We have also 
checked the superfluid density for the case that the magnon propagator V sharply peaked in momentum space at 
either q = or q ~ (t/o, ""/a, 7r/a). In each case it is again positive and it is somewhat smaller for the second case. 

The calculation of other response functions is similar to that for the Meissner effect and we give an example in 
Section VI. 



V. JOSEPHSON EFFECT BETWEEN ODD AND BCS SUPERCONDUCTORS 

Here we briefly consider the Josephson effect between odd and even (BCS) superconductors. To be specific, 
consider an S-I-S junction between an odd-frequency triplet superconductor on the left and an even-frequency triplet 
superconductor on the right. By applying the standard tunneling hamiltonian to this system, we find that the 
tunneling current has the form: 

</ ~ ^ \tk,kp\ 2 F odd(k, u)F even (kp, w) (5.1) 

k,kp,u; 

+ ^ l^,k' | 2 |^k",k'" \ 2 F odd (k, u)F odd (k", u)F even (kp, u)F even (k'" , w), 

k,kp,k",k"',o) 

where tk,kp is the hopping matrix element and we retain only the relevant fourth order term in the tunneling matrix 
element expansion. 18 For tunneling between either odd-odd or even-even gap superconductors the first term in the 
above equation is nonzero and dominates the Josephson current. The situation is drastically different in the present 
case. As long as there is no inelastic scattering in the junction, the leading term in the Josephson current between 
odd and even superconductor vanishes since the frequency integral in the first term is exactly zero.. However, the 
next term is nonzero and produces a current which corresponds to the tunneling of "pairs of pairs" with charge 4e. 
The very important implication of this result is for flux quantization. If we consider a combined ring, made out of 
even and odd superconductors, then the flux quantization inside this ring will be in units of 
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FIG. 9. Linearized gap equation in the presence of impurities. Solid lines are impurity averaged Green's functions. Dotted 
lines are electron-impurity interactions. 



The T c equation which replaces Eq. (3.14) is 



1 = VT 2 C Y, £( k ' w )£( k '' ^')^( k + k> + w')P(k + k', u + u') 



(7.1) 



k,k' 



where P(k + k', w + w') is the impurity vertex ladder seen in Fig. 9. It is given by 



P(k + k', w + w') = t[D (k + k') 2 + (w - u)')s] 



(7.2) 



Here, 1/r is the impurity scattering rate, D is the diffusion constant v^t/Z and s = [sgfn(w) — sgn(w')]/2. This result 
for the vertex is appropriate in the hydrodynamic regime |k + k'| < 1/vpT and w — w' < l/2r. Otherwise, P = 1. 

When Eq. (7.1) is evaluated with either of the magnon propagators mentioned in Section IIIC, i.e. the momentum- 
independent one or the spatially uniform one, finds a T c which is unchanged from the clean case up to corrections of 
order (1/E f t) 2 . 



Part of the work reported here was performed under the auspices of the Program on Strongly Correlated Electron 
Systems at the Los Alamos National Laboratory. The hospitality of the Aspen Center for Physics, where part of 
the work was carried out is acknowledged by EA and AVB. AVB was supported in part by a J.R. Oppenheimer 
Fellowship at Los Alamos. The support of NSF grants DMR92-21907 (EA), DMR92-25027 (DJS) and DMR92-22682 
(JRS) is also acknowledged. We are grateful to the following for helpful discussions and useful criticism: P. Coleman, 
A. Dugaev, CM. Eliashberg, L.P. Gor'kov, and P. Wolfle. 



1 A. Balatsky and E. Abrahams, Phys. Rev. B 45, 13125 (1992). 

2 P.W. Anderson, J. Phys. Chem. Solids, 52, 1313 (1991) 

3 F. Mila and E. Abrahams, Phys. Rev. Lett. 67, 2379 (1991). 

4 V.L. Berezinskii, Pis'ma Zh. Eksp. Teor. Fiz. 20, 628 (1974) [JETP Lett. 20, 287 (1974)]. 

5 E. Abrahams, A. Balatsky, J.R. Schrieffer, and P.B. Allen, Phys. Rev. B.47, 513 (1993). 

6 S. Yip and A. Leggett, unpublished 

7 A. Abrikosov, L.P. Gor'kov and I. E. Dzyaloshinskii, Methods of Quantum Field Theory in Statistical Physics, Dover, New 
York (1975). 

8 V.J. Emery and S. Kivelson, Phys. Rev. B26, 10812 (1992) 

9 A. Balatsky and J. Bonca, Phys. Rev. B48, 7445 (1993) 

10 P. Coleman, E. Miranda and A. Tsvelik, Phys. Rev. B49, 8955 (1994) 

11 Recall, W and A are related by A = W/Z , where Z is the self energy renormalization in the normal channel, which we take 
to be even in k and u> n . 

12 Note also that the odd gap cannot be obtained within the weak coupling BCS approximation which assumes no frequency 
dependence of the electron-electron interaction and therefore none in the gap A(k, u„). 

13 Should the first derivative vanish at r = 0, then it is necessary to go to the first non-vanishing odd r-derivative. In that case 
a much more complicated composite condensate containing more than one Cooper pair or more than one density fluctuation 
(or both) would enter. In our case, there is one of each. 

14 G.M. Eliashberg, Sov. Phys. JETP 11, 696 (1960). 

15 For example, J.R. Schrieffer, Theory of Superconductivity, (Benjamin, New York, 1984). 

16 J.R. Schrieffer, A.V. Balatsky, E. Abrahams, and D.J. Scalapino, J. Superconductivity, 7 501 (1994). 

17 J.R. Schrieffer and A.V. Balatsky, unpublished 

18 B.D. Josephson, Superconductivity, vol. 1, p. 423. Ed. R.D. Parks (M. Dekker, New York, 1969). 



VIII. ACKNOWLEDGEMENTS 



11 



